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We propose a simple experiment to study derealization and extinction in inhomogeneous bio- 
logical systems. The nonlinear steady state for, say, a bacteria colony living on and near a patch 
of nutrient or favorable illumination ("oasis") in the presence of a drift term ("wind") is computed. 
The bacteria, described by a simple generalization of the Fisher equation, diffuse, divide A — > A + A, 
die A — > 0, and annihilate A + A — > 0. At high wind velocities all bacteria are blown into an un- 
favorable region ("desert"), and the colony dies out. At low velocity a steady state concentration 
survives near the oasis. In between these two regimes there is a critical velocity at which bacte- 
ria first survive. If the "desert" supports a small nonzero population, this extinction transition is 
replaced by a derealization transition with increasing velocity. Predictions for the behavior as a 
function of wind velocity are made for one and two dimensions. 



Bacterial growth in a petri dish, the basic experiment of microbiology, is a familiar but interesting phenomenon. 
Depending on the nutrient and agar concentration, a variety of intriguing growth patterns have been observed [Q^Q. 
Some regimes can be modeled by diffusion limited aggregation, others by Eden models, and still others exhibit ring 
structures or a two-dimensional modulation in the bacterial density. At high nutrient concentration and low agar 
density, there is a large regime of simple growth of a circular patch (after point innoculation), described by a Fisher 
equation ||, and studied experimentally in Ref. Q|. 

Of course, most bacteria do not live in petri dishes, but rather in inhomogeneous environments characterized by, 
e.g., spatially varying growth rates and/or diffusion constants. Often, as in the soil after a rain storm (or in a sewage 
treatment plant), bacterial diffusion and growth are accompanied by convective drift in an aqueous medium through 
the disorder. By creating artificially modulated growth environments in petri dishes, one can begin to study how 
bacteria (and other species populations) grow in circumstances more typical of the real world. More generally, the 
challenges posed by combining inhomogeneous biological processes with various types of fluid flows || seem likely to 
attract considerable interest in the future. The easiest problem to study in the context of bacteria is to determine 
how fixed spatial inhomogeneities and convective flow affect the simple regime of Fisher equation growth mentioned 
above. 

A derealization transition in inhomogeneous biological systems has recently been proposed, focusing on a single 
species continuous growth model, in which the population disperses via diffusion and convection 0: the Fisher 
equation Q for the population number density c(x, t), generalized to account for convection and an inhomogeneous 
growth rate, reads j7||8) 




where D is the diffusion constant of the system, v is the spatially homogeneous convection ( "wind" ) velocity, and b is 
a phenomenological parameter responsible for the limiting of the concentration c(x, t) to some maximum saturation 
value (by competition processes of the kind A + A — ► Q). The growth rate U(jx) is a random function which 
describes a spatially random nutrient concentration, or, for photosynthetic bacteria, an inhomogeneous illumination 
pattern 0. If £/(x) is constant over the entire sample, then the convection term — v • Vc(x, t) has no effect on the 
growth of the bacteria. Only the introduction of a spatial dependence for the growth rate J7(x) makes the convection 
term interesting. In the following we consider the simple case of a "square- well potential" shape for J7(x), imposing 
a positive growth rate a on an illuminated patch ("oasis"), and a negative growth rate — ea outside ("desert") 
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I. INTRODUCTION AND RESULTS 




9c(x, t) jdt = L>V 2 c(x, t) - v ■ Vc(x, t) 
+f/(x)c(x, f) - 6c 2 (x,t) , 



(1) 




for |x| < -f , 
ea , for |x| > ^ , 



(2) 
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where W is the diameter of the oasis. Experimentally this could be realized using a very simple setup, which 
both illustrates the basic ideas of localization and derealization, and leads to interesting further questions. A one 
dimensional example is shown in figure where a solution with photosynthetic bacteria in a thin circular pipe, or 
annular petri dish, is illuminated by a fixed uniform light source through a mask, leading to a "square well" intensity 
distribution. The mask is moved at a small, constant velocity around the sample to simulate convective flow. (Moving 
the mask is equivalent to introducing convective flow in the system, up to a change of reference frame |L0|].) The 
bacteria are assumed to divide in the brightly illuminated area ("oasis") at a certain rate, but division ceases or 
proceeds at a greatly reduced rate in the darker region ("desert") outside. As a result, the growth rate in this 
continuum population dynamics model is positive in the oasis and small (positive or negative) in the surrounding 
desert region. Using this simple nonlinear growth model, we discuss predictions for the total number of bacteria 
expected to survive in the steady state, the shape of their distribution in space and other quantities, as a function of 
the "convection velocity" of the light source. 




FIG. 1. Experimental setup: a solution with photosynthetic bacteria in a circular pipe or a thin annular track in a petri dish, 
is illuminated only in a small area, while the rest of the sample is either kept dark or illuminated with reduced intensity. The 
light source is moved slowly around the sample to model convective flow. The bacteria are assumed to divide in the illuminated 
area ("oasis") at a certain growth rate a > 0, and die (or grow modestly) in the remaining area ("desert") with growth rate 
— ea. 

It is interesting to consider the class of biological situations discussed in this paper in the context of the " critical 
size problem" in population dynamics [jllj . In the critical size problem one asks for the minimal size of habitat for the 
survival of a population undergoing logistic growth and diffusion, where the region is surrounded by a totally hostile 
environment, i.e., no drift and an infinite death rate outside the oasis. We show here that the linearized version of the 
critical size problem is closely related to a well known problem in quantum mechanics, and present a generalization of 
this problem to include other types of surrounding environments, as well as the effect of drift. The linearized version 
of equation ([!]) around c(x, t) = reads 

dc(x.,t)/dt = £c(x,t), (3) 

with the linearized growth operator 

C = L>V 2 - v • V + [/(x) . (4) 

(We discuss later the validity of this linear approximation and compare the results with lattice simulations of the 
full nonlinear problem.) For nonzero convection velocity v, C is non-Hcrmitian, but it can still be diagonalized by 
a complete set of right and left eigenvectors, {<^(x)} and {0^(x)}, with eigenvalues r„ ]l2|,|7||, and orthogonality 
condition 
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d d x<^(x)<^(x) = ^ 



(5) 



(d is the dimension of the substrate, we focus here on d = 1 or d = 2). The time evolution of c(x,£) is then given by 

c(x, t)=^2 c n^(x) exp(r„i) , (6) 

n 

where the initial conditions and left eigenfunctions determine the coefficients {c n }, 

Cn = [ d d ^(x)c(x,i = 0). (7) 



Figure y shows the complex eigenvalue spectrum with the potential (g) for four different values of the convection 
velocity v, for a one dimensional lattice approximation to (|]) (see Appendix [A]) with periodic boundary conditions. 
The derivation of these results is discussed in section || below. At zero velocity C is Hermitian and all eigenvalues 
r„ are real. There are bound states (discrete spectrum) and extended or delocalized states (continuous spectrum). 
At finite velocities, all except one of the delocalized states acquire a complex eigenvalue. States with positive real 
part of the eigenvalue (ReT n > 0) grow exponentially with time, states with negative real part (ReT n < 0) decrease 
exponentially with time (see Eq. (||)). In a large one dimensional system the "mobility edge" fl3|| , which we define to 
be the eigenvalue of the fastest growing delocalized state, (i.e. the rightmost eigenvalue in the complex parabolas of 
figure ^J), is located at the overall average growth rate 

T* = (U)= { dxU(x)/L ~ -ea, (8) 



(up to corrections of order 0(1/ L) where the system size L is the mean circumference of the annulus in figure |l|). 
In figure |[ the eigenvalues of the localized states compose the discrete, real spectrum to the right of the mobility 
edge. With increasing velocity these localized eigenvalues move to the left by an amount proportional to v 2 /AD, and 
successively enter the continuous delocalized spectrum through the mobility edge, which remains fixed. The parabola 
broadens in the vertical direction - the imaginary parts of the eigenvalues of the delocalized states grow by an amount 
proportional to v. A given localized right cigenfunction (j)^ undergoes a "derealization transition" when the velocity 
reaches a corresponding critical derealization velocity v = u*, at which its eigenvalue r n has been shifted so far to 
the left that it just touches T* . At higher velocities it joins the parabola of eigenvalues describing a continuum of 
delocalized states. The ease with which such a derealization transition can be observed experimentally depends on 
whether there are growing delocalized eigenstates in the system, i.e. whether the mobility edge has a positive real 
value or not. 
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FIG. 2. Complex nonhermitian eigenvalue spectra (normalized by the difference of the growth rates inside and outside the 
oasis Uo = a + ea — 1) at velocities above and below the extinction transition. The spectra are extracted from numerical 
simulations of the lattice model described in Appendix [§], for a system of 1000 sites, oasis width W = 20 sites, D = 0.3, 
with growth rate U = —0.5 in the desert and U = +0.5 in the oasis, so that the average growth rate is is —0.48 (which is 
equal to the mobility edge F* up to finite size effects). The chosen velocity parameters (in units of the Fisher wave velocity 
in the oasis vf = 2a/ aD Q|) are v/vf ~ (circles), v/vf ~ 0.39 (plus), v/vf ~ 0.93 (triangles), and v/vf ~ 1.39 (squares). 
The point spectra are slightly offset in the y-direction so as to be able to distinguish the eigenvalues of the localized states 
for different velocities. As described in the introduction, the mobility edge remains roughly fixed, and the parabola of the 
delocalized eigenvalues opens up as v/vf is increased. The real, localized eigenvalues move to the left for higher velocities. At 
v/vf > 1.39 all states are delocalized. (This figure actually only shows the part of the spectrum which corresponds to the 
continuum problem. The lattice calculation also yields a left part of the spectrum - not shown here - which is an artifact of 
the discrete lattice.) 



In a large "deadly" desert (([/) ~ — ea < 0) all delocalized states die out, because the mobility edge lies to the left 
of the origin, as in figure |^. The growth rate of each localized eigenstate (fi^ then becomes negative at a corresponding 
"extinction" velocity v nc which is smaller than the corresponding derealization velocity w*. Thus, as convection 
is increased, the population dies out before it can delocalize. Later in this paper, we make specific predictions for 
the behavior of populations near the extinction transition, which occurs for v — vq c > v nc for all n > 0, when the 
eigenvalue of the localized "ground state" (fastest growing eigenfunction of £) passes through the origin. 

If the average growth rate (U) is positive {i.e. for a small enough desert or a small positive growth rate in an infinite 
desert), the mobility edge lies to the right of the origin and the derealization transition can indeed be observed at 
v — Vq where the "ground state" becomes delocalized. One expects to see universal behavior near this derealization 
transition, since there is a diverging correlation length in the system, which renders microscopic details irrelevant 
for certain quantities. We report predictions (see also iQ), for quantities such as the dependence of the localization 
length on the drift velocity as it approaches the derealization velocity, and the shape of the concentration profile 
near the transition. 

A special (universal) behavior is expected for the spatially average growth rate (U) = 0. In this case the der- 
ealization and extinction velocities coincide. Figure ^ summarizes the different scenarios in a sketch of the phase 
diagram for large systems with fixed well depth Uq = (a + ea), tuning the drift velocity, and the average growth 
rate (U) = —ea. Also shown in figure || is a horizontal transition line at (U) = separating a small velocity region 
(e < 0) where localized modes dominate the steady state bacterial population, from one (e > 0) containing a mixture 
of localized and extended states. The experimental signature of this interesting transition, (which could be accessed 
by increasing the light intensity for photosynthetic bacteria at fixed convection velocity) will be discussed in a future 
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publication jL4j . It is of course also possible to drive a population extinct at zero velocity simply by lowering the 
average growth rate. This special transition at (U) — —U c is indicated at the bottom of figure ||. 
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FIG. 3. Schematic phase diagram in one dimension for infinite system size, as a function of average growth rate (U) and 
convection velocity v. For a deep well (Uo = a + ea^> D/W 2 ) the extinction transition out of the localized phase occurs when 
v = for (U) = — U c , where U c — Uo- The diagram shows that if the growth rate is negative outside and inside the oasis 
(i.e. (U) < —U c ), then the only possible state is extinction at any velocity. If there is positive growth inside the oasis, but 
negative outside, a localized population can survive in the oasis, but only for small enough wind velocities v. Extended states 
are present for a small positive growth rate in the desert ({{/) > 0). In this case localized and delocalized states coexist for 
small velocities ("mixed phase"), while at large velocities all eigenstates are extended, as shown in figure []. The ground state 
becomes delocalized at the critical velocity v c which marks the phase boundary between the mixed and the extended phase. 

In section [n] we give details of the analysis of the one dimensional linearized problem for infinite and finite systems 
with periodic boundary conditions. In section [II some effects of the nonlinear term are discussed, especially for 
experiments near the extinction transition, and in section [i^ the two dimensional case is discussed. The appendices 
contain some details on the analytic computation of finite size effects (Appendix ^), a brief discussion of a lattice 
model |Q corresponding to the analytic continuum theory ( Apppendix |b|) , and a discussion of dimensionless quantities 
measurable in experiments (Appendix ^) . 



II. LINEARIZED GROWTH IN ONE DIMENSION 



If the left and right eigenfunctions (f>^' L (x.) are localized (i.e. if the convection velocity is small enough, so that 
4>n' L ( x ) decays exponentially with the distance from the oasis) , one may eliminate the convective term in Eq. (Q) via 
the transformation 

<^ L (x) = exp(±v • x/2Z>)V„(x) (9) 

(+ refers to the right eigenvectors and — to the left eigenvectors). The eigenvalue equation associated with the 
linearized growth operator (||) becomes Hermitian |1| 

r„V«(x) = BV 2 i(x) + C/(x)V„(x) - (v 2 /4Z>)W,(x) , (10) 

and is equivalent to the familiar square well potential problem much studied in quantum mechanics. With the 
identifications a + ea = Uq/Ji, where Uq is a quantum well depth, and fi is Planck's constant, T n + ea + v 2 /4D = \E\/h, 
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where E is a quantum energy level, and D = h/2m, where m is a mass in the equivalent quantum problem, we can use 
well known quantum mechanical results |l6|]l7f l . The left and right eigenfunctions at finite velocity are then related to 
the eigenstates of the Hermitian problem (|10|) via the transformation (||), while the eigenvalues undergo a rigid shift 

r„(«) = r n (v = o) - (ii) 



A. An Oasis in an Infinite Desert: Localized Populations and the Extinction Transition 

In an infinite one dimensional system, localized solutions for </>^ ,L (x) are given by Eq. (0) with [Q 



' Ai 7n exp(K„x) , for x < — ^ 



ip„(x) = < 



2 ' 



e^qy(ik n x) + £?2,n exp(— ik n x) , for — — < x < — , (12) 



k A 2 , n exp(-K„x) , for x > 2 , 



where vl2,n, B\ n , and -B2,n are constant coefficients, and 



Kn = y/(T n (v = 0) + ea)/D, (13) 

and 



K = v /(a-r n (u = 0))/£>, (14) 

as can be seen by substituting the above Ansatz for ip(x) into Eq. (|To|). Eq. (|l3| ) implies 

r„(w = 0) = L>k^ - ea. (15) 

To compute T n (v = 0) one matches both ip n (x) and dt/j n (x) / dx at X = ±W/2, which determines the coefficients 
2, up to an overall multiplicative factor, as well as the eigenvalues {T n (v = 0)}. When solving the 
Hermitian problem one may use the fact that all the eigenfunctions admit a well defined parity, i.e., they are odd or 
even under the transformation x — > —x. Even integers n correspond to bound eigenstates with even parity, where 
if) n (x), is symmetric under x — > —x. In such a case one obtains the eigenvalue equation for the quantity 

C„(r„) = k n W/2 = W/2^(a - T n (v = 0))/D , (16) 

namely, 



COt(Cn) = CnX/y/1 ~ (CnX) 2 , (17) 
(which is equivalent to cot(k n W/2) — k n /n n ), with 



x = 2^D/{a + ea)/W = 2/{^k 2 n + n 2 n W) . (18) 

The dimensionless parameter x, measures the ratio of kinetic to potential energy in the equivalent quantum problem. 
For bound eigenstates with odd parity ip(x), (^>(x) antisymmetric under x — > —x, denoted by odd n ) one obtains 



COt(C„) = -Vl-(CnX) 2 /(Cn*) ■ (19) 

The highest eigenvalue (with its nodeless, positive eigenfunction ip (x)), corresponds to the largest growth rate i?e{r }, 
and is therefore expected to dominate the system in most cases at long times, as seen from Eq. (^). Its eigenvalue 
condition ( pT| ) for n = leads to fill 

r o (w = 0) + ea= (a + ea)f(x), (20) 

where f(x) is a monotonically decreasing function such that f(x) ~ 1 — 7r 2 i 2 /4 for x <C 1, and f(x) ~ 1/x 2 for x ^> 1. 
Upon inserting Eq. (p0|) into the expressions for k„ and k n one obtains 
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Ko = ^{a + ea)f{x)/D (21) 

and 



fc = y/((a + ea)0—f(x)))/D . (22) 

When i < 1, the "potential well" is very deep, and one finds the usual particle in a box result for ro(u), with 
correction term proportional to v 2 arising from the change of variables <M) and (EjJ), namely 



T a (v) ~ (v 2 - v 2 )/(AD) =a- Dn 2 /W 2 - v 2 /4D , (23) 

with 



v c = 2Dy/(a/D - tt 2 /W 2 ) , (24) 



k = ^{a + ea)/D-TT 2 /W 2 (25) 

and 

k a =7r/W. (26) 

For x 1 one finds 

r (w) (v 2 - w 2 )/(4 J D) = (a + ea) 2 (W/2) 2 /D - eo - D(v/2D) 2 (27) 

with 



w c = 2iV((a + ea)W/2) 2 /-D 2 = ea/- , (28) 



Ko = V(a + ea) 2 (Wy2) 2 AD 2 (29) 

and 



fco = y/(a + ea)(l - l/x 2 )/D. (30) 

If we take as an effective diffusion constant for motile bacteria D = 6 • 10 _6 cm 2 /sec, and a growth rate a — 10~ 3 
/sec in the oasis and a much smaller growth rate outside (0 < |e| << 1), we get for a W = 2cm diameter oasis, 
x = 0.077 -C 1, and an "extinction velocity" v c ~ 1.5/x/sec, which is comparable to the Fisher wave velocity || in the 
oasis, vp = 2\J aD = 1.5/x/sec. 

B. Finite Size Effects and the Delocalization Transition 

An experimental finite system with periodic boundary conditions is depicted in figure [l]. The eigenvalue equation 
for fc n (r„) and fc n (r n ) of the corresponding linearized problem of an oasis of width W in a finite desert of extent 
L > W (with L being the mean circumference of the circular region in figure [j]), is obtained using the Ansatz in 
Eq. (§) with 

(Ai t n exp(K„x) + A 2 ,„ exp(-K„x)) , for ~ < x < , 

ip n {x) = { (B^ n exp(ik n x) + B 2 , n exp(-iknx)) , for -\ < x < \ , (31) 

(Cx, n exp(K n x) + C 2 , n exp(—K n x)) , for \ < x < j , 

and matching (j)^' L (x) and d(j)^ ,L (x)/dx at the edges of the well and at the edges of the sample (imposing periodic 
boundary conditions). One finds the eigenvalue equation Q (see also Q) 

2 k„n n (cosh(Lv/(2D)) - cos(k„W)cosh(K n (L - W))) 

+ (fc 2 - Kl)sin(k n W)sinh(K n {L - W)) = . (32) 
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For v/(2D) < Re{n n }, and large L, equation ( |3^ ) yields 

cxp(L(v/(2D) - n n )) = (33) 
cos(fc„VF) - k n /(2n n ) ■ (1 - {n n /k n ) 2 ) sm{k n W) . 

The left hand side vanishes in the limit L — > oo, and the equation reduces to the bound state equations of a single 
square well (Eq. ( fL7| ) for even parity solutions and Eq. ( [l9| ) for odd solutions). For finite L, and v/2D < Re{K n }, the 
deviation of the "localized" solutions n n and k n from their u = values for small v is exponentially small in L p2| . 
These "localized" or "bound state" solutions, are characterized by an exponential decay of the bacterial density in 
the desert with a correlation length 

£„ ~ (Re{ Kn } - v/(2D))- 1 (34) 

and T n (v) strictly real. However, "delocalized" or "scattering" solutions also exist, with T n (v) complex, and nontrivial 
dependence of n n and k n on v, even in the limit of large L. As the velocity is increased, the nth localized eigenstate 
becomes delocalized (£ n — * oo) at the critical derealization velocity given in an infinite system by 

v* = 2DRe{n n } . (35) 

This implies 

~ l/« - vY (36) 

with the (universal) critical exponent v = 1. We saw that with increasing velocity, the growth rate Re{T n (v)} for a 
given eigenstate decreases. It becomes negative above the corresponding extinction velocity v nc . We therefore expect 
that the derealization transition for the "ground state" (which tends to dominate the long time behavior) can be 
observed only if Vq < vq c . In the following we discuss the three desert scenarios, e > 0, e < 0, and e = 0. 

(1) For e > (a "deadly" desert), of big enough size L, one finds that all delocalized states die out exponentially 
with time (v nc < u*). The population is localized around the oasis at small v and extinct at high v. Figure ^ is a 
plot of the eigenvalues T(v) in the complex plane for this case, as derived for the lattice model discussed in Appendix 
[B|. The lower part of figure || shows a series of profiles of the ground state eigenfunction close to the extinction 
transition. In small enough systems, such that the total effective growth rate f Q L U (x)dx is positive, delocalized states 
can actually have a positive growth rate even for e > 0. (See also Eq. (|A2[ ) in Appendix [a|, with 6k ~ 0(1/ L), and 
R ~ 0(1/ L) .) In this case, the system is small enough so that the bacteria can traverse the desert quickly, and on 
average won't die before reenterring the oasis in a circular pipe. 

(2) If e < 0, delocalized states should be observable even for very large systems, because the "desert" can support 
modest growth, although at a much smaller rate than in the oasis if |e| <C 1. Growing delocalized eigenstates are 
present, even for v = 0, and the population is a superposition of fast growing localized states and more slowly growing 
delocalized ones. As a drift velocity v increases, the n'th localized eigenstate delocalizes at v = u* with a positive 
growth rate ReT n (v^ ) (i.e. w* < v nc ). This case allows for an experimental observation of the derealization transition: 
as the velocity is increased, more and more eigenstates delocalize. The eigenvalue spectrum for two different values of 
v is shown in figure [| (a) and (b) . We can see that the spectrum at the derealization transition is slightly different 
depending on whether an "even" eigenfunction or an "odd" e igenf unction is about to delocalize next ("even" and 



"odd" are to be understood in the sense explained in section II A ). If an odd eigenfunction is about to delocalize 
next, there exists a delocalized state which has a purely real growth rate (at the tip of the parabola in the spectrum 
of figure p(b)), while no such state exists when an even eigenfunction is about to delocalize, as in figure |(a). The 
essential characteristics of the spectrum are derived in Appendix 

For experiments, we focus on the derealization of the ground state, since it is the fastest growing eigenstate, which 
dominates the system near the oasis. The ground state delocalizes at the highest derealization velocity v = Vq, 
i.e. for v > Vq all states are delocalized. The lower part of figure |^ shows a series of profiles of the ground state 
eigenfunction close to derealization. One sees that at the derealization velocity Vq, the correlation length £o reaches 
the system size. In an infinite system it diverges as in Eq. ([36]). Near the derealization transition, we expect certain 
quantities to become independent of microscopic details of the system. These universal quantities will only depend on 
general properties, such as symmetries, dimensions etc.. Simple models that only share these general properties with 
the experimental system will be sufficient for precise predictions for these quantities near the transition. An example 
of such a quantity is the critical exponent v in Eq. (|36|) . Details for the square well system and more general random 
systems will be presented in a future paper jljj . 

(3) For e = 0, the growth rate of the bacteria exactly balances the death rate in the desert, they only diffuse. In this 
case the derealization velocity and the extinction velocities coincide: w* = v nc in an infinite system. There is again 
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a diverging correlation length leading to universal critical behavior near the derealization transition |l4|] . In a finite 
system the average growth rate dxlI(x)/L is positive. In small systems, one therefore expects to see delocalized 
states with a positive growth rate in experiments for this case as well. 

These results for one dimension in the limit L — ► oo are summarized in figure ^ which shows a schematic phase 
diagram for fixed well depth Uo = (a + ea), as a function of the convection velocity and the average growth rate 
(U) = J dxlI(x)/L, in the limit L — > oo. 
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FIG. 4. A series of steady state population profiles for a system with negative average growth rate, at different drift velocities, 
for the linear and the nonlinear case. The x coordinate is normalized by the width W of the oasis. In the nonlinear case the 
populations c(x,v) are divided by the population at zero velocity in the middle of the oasis c(0, 0). In the "linear" case, we 
simply assume the validity of Eq. (0) and rescale the curves such that their maximum values match those of the nonlinear case. 
The profiles were extracted from a lattice model with 1000 sites, an oasis of width W — 200 sites, with growth rate U = 0.5 
inside the oasis and U — —0.5 outside, and diffusion constant D = 30. The velocity parameters are taken to be v/vf — (solid 
line), v/vf = 0.93 (dotted), v/vf = 0.98 (dashed), and v/vf = 0.99 (dash-dotted). Close to the extinction transition, where 
the population is small, the agreement between the nonlinear and the linear solution becomes quite good. The oasis occupies 
the region \x/W\ < 0.5 in the figures. 
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FIG. 5. Sketches of complex nonhermitian eigenvalue spectra at two different velocities: (a) velocity v a with v? < v a < v£, 
where n = 6 denotes an even eigenstate, i. e. all eigenstates with n > 6 are delocalized and those with n < 6 are localized; and 
(b) at a higher velocity v b > v a with v£ < v b < vl where n = 5 denotes an odd eigenstate. Upon increasing the velocity from 
v a to Vb the n = 6 eigenstate becomes delocalized at v = vl . 
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Bacterial Density near Derealization 

Nonlinear Case 




x/W 

FIG. 6. A series of steady state population profiles for a system with positive average growth rate, at different drift velocities, 
for the linear and the nonlinear case. The x coordinate is normalized by the width W of the oasis (|a;/VK| < 0.5). In the nonlinear 
case the populations c(x,v) are divided by c(0,0), the population in the middle of the oasis, at zero velocity. Curves for the 
"linear" case, were computed as in figure m The profiles were extracted from a lattice model with 1000 sites, an oasis of width 
200 sites, with growth rate U = 0.1 inside the oasis and U = 1.1 outside, and large diffusion constant D = 300. The velocity 
parameters v/vf are taken to be v/vf = (solid line), 0.17 (dotted), 0.33 (dashed), 0.50 (long-dashed), and 0.83 (dash-dotted). 
In both cases the population becomes delocalized at high velocities with increasing velocity when the correlation length reaches 
the size of the system. 

C. Delta- Function Potential Well 

Results for a delta-function-like oasis in one dimension, can be easily derived as a special case of the square well 
potential, by taking the limit e — > 0, a — > oo and W — > with aW = const = Vq. In this case Eq. (^) simplifies to 

n(cosh(Lv/(2D)) - cosIi(kL)) + V smh(KL)/D = , (37) 

which is the same equation as obtained in jl2j. With the identification k = {—i)K the results derived there can be 
applied here: the delocalization picture is the same as in figure §(a), except that there exists only one (even parity) 
bound state solution for i?e{«} > v/(2D). Again two critical velocities v c q and Vq emerge, in accordance to the above 
discussion for a square well potential. 

III. EFFECTS OF THE NONLINEARITY 

We can also estimate the effects due to the nonlinear term in equation (|l|), which leads to a saturation of c(x, t) for 
r„(u) > 0. The equation of motion becomes 

dc(x, t) jdt = £c(x, t) - bc 2 {x, t) , (38) 

with C given by Eq. (^). (The coefficient b can be set to 1 by rescaling the density c(x,t) by 6, as in Appendix [c]. ) 
The solution can be expressed in terms of the complete set of right eigenstates for the linear problem with new time 
dependent coefficients c n (t) with 
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c(x,t) =^4)^(it) (39) 

n 

and 

dc n (t)/dt = T n c n (t) - ^ 

tV n mm ,' C m (t)c m ,(t), (40) 

m,m' 

where the mode couplings are given by [Q 

w n , mm > =b dx<t>%(x)<j>*(x)<j>*,(x) . (41) 



In general one expects that through the mode couplings the fastest growing eigenstate suppresses the growth rate of 
the other eigenstates, provided the corresponding couplings are large enough J^JlJl- In the mixed phase of figure || we 
expect that the fastest growing bound state suppresses the other bound states, and the fastest growing delocalized 
state suppresses the other delocalized states H . 



A. Effects of the Nonlinearity at the Extinction Transition 

For e > (and large enough systems so that J Q L dxU(x) < 0), with velocity v just below the extinction velocity 
t>0c, we expect that the growing ground state term co{t)4>o( x ) dominates the summation (^9|) for c{x, t). In a first 
approximation we neglect all c m (t) with m > and find 

dc (t)/dt = r oCo (t) - w c$(t), (42) 

with wq = wo, oo = bj d d x(j)Q (x)(4*q(x)) 2 > 0. At long times 

^^' H.M^M-.Wr. (43) 

with asymptotic behavior lim^oo Co(t) = Tq/wq. Thus, the steady state population profile should be given approxi- 
mately by 

c*(x) =T <l > X(x)/wo, (44) 

where 

c*(x) = limt->aoc(x, t) . (45) 
The total steady state bacterial population Nq is 

«b (46) 

It follows that 

N ~ (r /tflo) / dx^\x) . (47) 



Since To ~ (v — v c )(v + v c ), and J dx<j>Q'(x) ~ const + 0(v — v c ) as v — * v c from below, one finds that near the 
extinction transition 

{v c — v , for v —>■ v~ , 
(48) 
, for v > v c , 

Figure ^ shows the total population Nq as a function of the convection velocity v for an oasis in a desert with negative 
growth rate, which is large enough so that the average growth rate is negative also. The data was obtained from a 
lattice model for the nonlinear problem (see Appendix [A]) . The displayed convection velocities range from v — to 
velocities larger than the extinction velocity v > v c . One can see that Nq decreases linearly with v — v c for v — * v~ , 
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as predicted in Eq. (|4q). Figure |] shows a series of population profiles for increasing velocity, for the linear case (i.e. 
taking the ground state of the linear problem as a solution for the steady state), and for the nonlinear case. One can 
see that close to the extinction transition the ground state of the linearized growth is an excellent approximation, as 
is to be expected since there the bacterial density c(x) becomes small, so that mode couplings other than wq induced 
by the nonlinear term be 2 (x) become small relative to the linear terms. 

Total population versus drift velocity 



(Extinction transition) 




FIG. 7. Numerical results for the total number of bacteria No(v) (normalized by the total number at zero velocity No(0)) as 
a function of v/vp, where vf = 1\faD is the Fisher wave velocity. The data were obtained using a 1 dimensional lattice model 
with 100 grid points (see Appendix [§]) . The diffusion constant is D — 3., the width of the oasis is W = 20, the growth rate 
is -0.9 outside the oasis and 0.1 inside. The coefficient of the nonlinear term in the equation of motion is b = 1. The average 
growth rate is —0.7 < for this system, so that vo c < t>o, and the system goes through the extinction transition at v/vf — 0.65, 
above which Nq(v) vanishes. The approximation v c ~ vf mentioned in the text does not work quite as well here as in other 
cases because the condition i < 1 necessary for this approximation is violated weakly for the model parameters chosen here. 
(x ~ 0.35). 

B. Limit v — > 

In the linearized case, because Eq. ( pToj ) is even in v, we expect No(t) = J d d xc(x, t) ~ const — const' ■ v 2 + 0(v 4 ) for 
small velocities. In the nonlinear case this symmetry is broken, because of the v dependence of the coefficients w nm ,m' in 
Eq. (fill ) which arises from the transformation (^). One then expects in the steady state No ~ const—const' ■ |v| +0(v 2 ) 
for small velocities. The constants depend on U(x) and other nonuniversal details. 

C. Effects of the Nonlinearity at the Derealization Transition 

Figure |^ shows the same information as figure || for the derealization transition for an oasis in a background with 
positive growth rate. One can see that the linear approximation (i.e. taking the ground state of the linear problem as 
a solution for the steady state), becomes best for high velocities [Q, when the drift term in the equation of motion 
becomes dominant compared to the nonlinear term. Furthermore the nonlinear solution at low velocities is in the 
"mixed phase", i.e., it is a superposition of extended and localized eigenstates of the system. A decomposition of 
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the v = nonlinear solution into the eigenstates of the system shows |14J] that its leading contributions are from 
the localized ground state eigenfunction and the fastest growing delocalized eigenfunction. As the drift velocity is 
increased, the contribution of localized eigenstates clearly diminishes. Above the derealization transition, the steady 
state is composed only of delocalized modes. 



IV. TWO DIMENSIONS, INFINITE SYSTEM SIZE 

Much of the above analysis can be adapted to two dimensions. We again use the transformation (^|) and consider 
a circular oasis of diameter W in Eq. (|l0|). The analysis is straightforward, so we omit most of the details. Following 
reference ]l2| ], the qualitative behavior of non-Hermitian eigenvalue spectra for linearized growth in two dimensions 
should be as follows: As in one dimension, localized states associated with the potential well will lie in a point spectrum 
on the real axis. Extended states, however, will occupy a dense region in the complex plane with a parabolic boundary 
(see figure ||). With increasing w, the localized point spectrum will again migrate into the continuum. 
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FIG. 8. Sketch of a complex nonhermitian eigenvalue spectrum in two dimensions with localized and delocalized eigenstates. 
The sketch uses a computation of the spectrum for a delta function potential in two dimensions J^] , in which case there is only 
one bound state (to the right of the imaginary axis). We expect qualitatively the same result for the finite two dimensional 
well discussed in the text, except that multiple bound states with real eigenvalues will occur, just as in one dimension. 

In the localized regime the spatial distribution of bacteria is in the linear approximation given by the convection- 
distorted ground state eigenfunction, 



V ( x ) 



d i oe X p(±vx/2£>)J (fc n |x|) 
C 2 ,oexp(±vx/2L>)iY ( 

K n |x| 



for |x| < W/2 . 
for |x| > W/2 , 



(49) 



where Jq(x) and Kq(x) are Bessel functions of order zero [E0|. The constants Ci^o and 6*2,0 are chosen such that 



L R.L 



(x) and its radial derivative are continuous at |x| = W/2. The lowest energy eigenvalue is again given by eq. (pC 



but with different results for f(x) |16|: for x <C 1 one finds f(x) ~ 1 



for x 3> 1 one finds f{x) ~ x 2 exp(— Ax 2 ). The behavior of the ground state eigenvalue is thus 

r ~ a - AD(2A05) 2 /W 2 - Dv 2 /{AD) 

for i<l and 



where a\ ~ 2.405 is the first zero of Jo(a); 



(50) 
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r ~ D{A/W 2 ) exp[-16£>/(VF 2 (a + eo))] - D(v/2D) 2 



(51) 



for x 1. From these results one can compute vq c for the 2 dimensional case. As in one dimension, one again finds 
for e > ("deadly" oasis), and large enough systems, that v c q < uj-j, where Vq/2D ~ i?e{«o}, and v$ c is the velocity 
above which the ground state growth rate Tq becomes negative. 

Figure ^| shows a contour plot for the spatial distribution of bacteria in the linearized case at convection velocity 
v below Vq c in 2 dimensions, obtained from equation (|l9|). Although the full nonlinear two dimensional problem 
seems tractable numerically, we expect that close to the extinction transition the curves for the linearized case will 
approximate the shape of the steady state, because the ground state dominates sums like that in equation (B9h . 




FIG. 9. Contour plot of the distribution of bacteria in two dimensions as obtained for the linearized problem for convection 
velocity v = 0.2fj,/sec < v c = 1.5/i/sec (directed to the right), with D = 6 • 10 _6 cm 2 /sec, and growth rate a — 10 _3 /sec 
in the oasis (of diameter W = 2cm) and a much smaller growth rate outside (|e| = 0.0001 1). The dashed line indicates 
the circumference of the oasis, the coordinates are given in units of the radius of the oasis. At high velocities the linear 
approximation is expected to give the steady state generated by the full nonlinear equation, provided that the ground state 
dominates the long time behavior in Eq. (E9I) (see text). 



APPENDIX A: COMPLEX NONHERMITIAN EIGENVALUE SPECTRA FOR FINITE SYSTEM SIZE 



To compute t he g rowth rate T n (v) near the derealization of the nth eigenstate, we set k = v/(2D) + Sk~ in where 
liniL^oo 8k = ]12fl, and take v to be close to the corresponding v*. Upon expanding the right hand side of equation 



15 



( |34| ) in powers of v — u* and 6k — in, one obtains: 

exp(— L(Sk — ik)) — c n {v — v^)/v + 0(6k — in) 

+ 0((v-v* n f), (Al) 

where c„ is the derivative of the right hand side with respect to v at v = v * and 5k— in = 0. For even n one finds c„ < 0, 
and for odd n one finds c„ > 0. Hence, for even n, and v < u*, one has to leading order 5k ~ ln[v / (c n (v — v*))]/L 
and k = 7r(2m + l)/£, with to any integer. For v > u*, 5k ~ ln[w/(c„(?;* — v))]/L and re = 2-Km/ L. Similarly, for odd 
n and v < v*, 5k ~ ln[w/(c„(w* — and re = 2irm/L, and, for u > u*, one has 8k ~ ln[w/(c„(u — w*))]/L and 

k = 7r(2m + At v = w* in either case, 5k ~ 0((lnL)/L) and R ~ 0{{Yn L) / L) ]T^| . For given u and k the value 

of r results from above definition of «, r(u) = Z)k 2 — ea — D(v/(2D)) 2 , and thus 

I» = Z?(5k + w/2Z?) 2 - Hk 2 - D(v/(2D)f - ea 

-2i£>(oK + v/2L>)k. (A2) 

These results are illustrated in figure 0. 



APPENDIX B: NUMERICAL ANALYSIS OF A DISCRETE LATTICE MODEL 



A discrete lattice model, originally inspired by the physics of vortex lines has proven very helpful for the 
numerical analysis of our problem. The corresponding lattice discretization (with lattice constant Zo) of the nonlinear 
equation in d dimensions reads M 



dc x (t) w 



^[e g - e "c x+e „(<) + e- g ' e "c x _ e „(t) - 2cosh(g ■ e,)c x (t)] 



dt 2 

+?7(x)c x («) - k£(t) , (1.18) 



i/=i 

2 / 



where c x (t) is the species population at the sites {x} of a hypercubic lattice, and the {e^} are unit lattice vectors. 
Furthermore, (w ~ D/£q), where D is the diffusion constant of the corresponding continuum model, and g ~ v£q/(2D), 
where v is the convective flow rate of the continuum model. U (x) and 6 have the same interpretation as in the 
continuum model @. The subtraction in the first term insures that c x (t) is conserved (4^ J2 X c x(t) = 0) if U(x) = b = 
0. There are two constraints on the mesh size or lattice spacing £q, which ensure that the model correctly describes 
the continuum limit for a given eigenstate <f>^, derived from the condition that £o|V<^(x)|/<^„ (x) <C 1. For small v 
one requires 



k n £ < 1 and K n £ < 1 . (Bl) 



For high velocities the condition becomes 



v£ /(2D) « 1 , (B2) 
as follows from Eq. (pM. The lattice simulations shown in this paper are well within these limits. 



APPENDIX C: DIMENSIONLESS QUANTITIES 

The equation of motion can be rewritten in dimensionless form, by introducing rescaled coordinates, y = x/W, 
where W is the width W of the oasis, and rescaled population densities c = c/c s , where c s is roughly the saturation 
value of the bacterial density in the oasis (up to 0(e)): c s — a/b (see Eq. ([]])). One then obtains 

(W 2 /D)dc(x, t)/dt = V 2 c(x, t)-v- V y c(x, t) 

+4[(C/(x)/f/ )c(x, t) - c 2 (x, *)] , (CI) 

where v = vW/D is the dimensionless rescaled drift velocity, and vf = 2yJW 2 a/ D is the dimensionless rescaled 
Fisher wave velocity vf — 2\J Da in the oasis, which gives the speed at which a Fisher wave |^| would propagate in 
the oasis. The velocity vp also gives a rough estimate for the velocity at which the extinction transition takes place, 
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in the appropriate parameter regime (e < 0) of the phase diagram of figure [|. The basic time scale of the system is 
set by the diffusion time W 2 /D, which is the time it takes a bacterium to diffuse across the oasis. 
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